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Abstract 

The classical relativistic wave equations are presented as partial difference equa- 
tions in the arena of covariant discrete phase space. These equations are also 
expressed as difference-differential equations in discrete phase space and contin- 
uous time. The relativistic invariance and covariance of the equations in both 
versions are established. The partial difference and difference-differential equa- 
tions are derived as the Euler-Lagrange equations from the variational principle. 
The difference and difference-differential conservation equations are derived. Fi- 
nally, the total momentum, energy, and charge of the relativistic classical fields 
satisfying difference- differential equations are computed. 

PACS: ll.lO.Ef; ll.lO.Qr; ll-15.Ha 
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1. Introduction 

Partial difference equations have been studied [1] for a long time to investigate 
problems in mathematical physics. Moreover, modern numerical analysis 
[2], which studies the differential equations arising out of various physical 
problems approximately, is based upon finite difference (ordinary or partial) 
equations. 

Recently [3], we have formulated the wave mechanics in an exact fash- 
ion in the arena of discrete phase space in terms of the partial difference 
equations. This new formulation includes (free) classical relativistic Klein- 
Gordon, Dirac, and gauge field equations. 

The proofs of the relativistic invariance or covariance of various partial 
difference and difference-differential equations are quite subtle. However, we 
have managed to provide such proofs in Section 4. 
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The Euler- Lagrange equations for the partial difference and difference- 
differential equations and Noether's theorems for difference and difference- 
differential conservation laws are derived in Appendix I and Appendix II 
respectively. 

Finally, the total momentum, energy and charge are computed for vari- 
ous relativistic classical fields obeying difference- differential equations in Sec- 
tion 5. (We have not computed conserved quantities for wave fields satisfying 
partial difference equations for a physical reason to be explained later.) 

Moreover, the stress-energy-momentum tensor and the consequent total 
momentum-energy as computed from a general Lagrangian are somewhat 
incomplete in this paper. However, in the next paper, these quantities for 
the free Klein-Gordon, electro-magnetic, and Dirac fields are furnished with 
exact equations. 



2. Notations and preliminary definitions 

There exists a characteristic length I (which may be the Planck length) in this 
theory. We choose the absolute units such that h = c = i = 1. All physical 
quantities are expressed as dimensionless numbers. Greek indices take from 
{1,2,3,4}, roman indices take from {1,2,3}, and the capital roman take 
from {1,2}. Einstein's summation convention is followed in all three cases. 
We denote the flat space-time metric by rj^u and the diagonal matrix [rj/^u] '■— 
diag[l, 1, 1 — 1]. (The signature of the metric is obviously +2.) We denote 
the set of all non-negative integers by N := {0} U {Z^} = {0, 1, 2, 3, . . .}. An 
element neNxNxNxN and an element (n, t) e x M can be expressed 
as 

n = (r^^r^^r^^r^^), n" eN, // e {1, 2, 3, 4} ; (lA) 
(n,t) = (n^n^n^^), eN, jG {1,2,3}, teR. (IB) 

Here and elsewhere, a bold roman letter indicates a three-dimensional vector. 
In this paper, the equations in the relativistic phase space are denoted by 
(..A), whereas the equations in the discrete phase space and continuous time 
are labelled by (..B). Both formulations are presented up to the difference and 
difference-differential conservation laws. Subsequently, only the difference- 
differential equations are pursued. The physical meanings of the quantum 
numbers are understood from the equations 

ixn' + ilf)'^2n'^ + l, /.e {1,2, 3, 4}, (2) 

where x'^ denote the space-time coordinates and stand for four-momentum 
components in quantum mechanics. 
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FIG. 1 One discrete phase plane. 



Therefore, gives rise to a closed phase space loop of radius \j2n^ + 1 in 
the /x-th phase plane. Field quanta reside in such phase space loops where the 
measurements of angles are completely uncertain (see Fig. 1). Phase space 
loops can be interpreted as degenerate phase cells. 

We shall encounter three dimensional improper integrals in computing 
the total conserved quantities. Those integrals are always defined to be the 
Cauchy-principal- value: 



/(k)rf3k:= lim 

L— +0O 



L L L 



f{ki, k2, ks) dkidk2dk3 



-L -L -L 



(3) 



Let a function be defined by / : — ^ IR (or / : N'' ^ C). Then the 
right partial difference, the left partial difference, and the weighted-mean 
difference are defined respectively by [3,4]: 



A^/(n) := /(..,n'^ + !,..)-/(..,<..), 
A;/(n) :=/(..,n^..)-/(..,n^-l,..). 



(4i) 
(4n) 



A*/(n) := (1/V2) [V^^ij^ f{..,n^ + 1, ..) - fi..,n^ - 1, ..)] (4iii) 
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It is clear that the partial difference operators A^, AJ|^, are all linear and 
the operators — iA* are self-adjoint [3]. By direct computations, we can 
prove the following generalizations of the Leibnitz rule: 



AMn)9{n)] = /(.., + 1, ..)A,g{n) + g{n)A,f{n) , 



(5i) 



K[fi^)9in)] = f{n)A'^g{n) + g{.., n>^ - 1, ..)A'J{n) , (5ii) 



A*[/H5(n)] 



f{..,n^^+l,..)A*g{n)+g{..,n^-l,..)A*f{n) 
-f{..,n'^ + l,..)g{..,n'^-l,..)A*{l), 

/(.., ..)A*g{n) + g{.., + 1, ..) A*/(n) 

-/(.., n^-l,..)y(..,n'^+l,..)A*(l), 



(5iii) 



A*(l) := (1/^2) fV^^^^TI- Vn^l , limA*(l) = 0, (5iv) 



n**— >oo 



A/, { [0(n)V^(.., n^' - 1, ..) + n'^ - 1, ■OV'N]} 
^^/2[(t>{n)A*^l;{n)+i;{n)A*(t>{n)]. 



(5v) 



In the left-hand side of the equation (5v), the index fi is not summed. 

We shall furnish a few more rules involving finite difference operations in 
the following equations: 



A,(t){n) + V^A'Xn) = \/2 [A*0(n) - 0(n) • A*(l)] , (6i) 



^/^^FTl A,A*(t>{n) + ^/^A',A#(/)(n) 
= v^[A*A*0(n)-A*0(n)-A*(l)], 



(6n) 



V^^FTilA^m? -V^[A'^4^{n)f 

= ^{A*[0(n)]2 - 20(n)A*0(n) + [^(n)]^ • A*(l)} , 



(6in) 
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= ^/2 {A*(A*0 ■ AJ0) - A*0 ■ (A*A#0) - A*0 ■ (A*A^ 
+(A*0)-(A*0)A*(1)}, 



(6iv) 



V^^^{A^A,){A^B,) -V^{A'^A,){A'^B,) 

= V2 [A*{A,B,) - A,A*B, - {A*A,)B, + A,B,A*{1)] 



(6v) 



Here, neither the index /j, nor the index v is summed. 
Now the rules for the summations will be listed. 



Y: A,f{n) = /(.., + 1, ..) - /(.., TVf , ..) , 



(7i) 



ATM 



E (A;/(n) = /(..,iV2^..)-/(..,iVr-l,..), 



(7ii) 



^2 r 

E A*/(n) = (1/V2) 0Vrri/(..,Ar|' + l,..) 



-ViVr/(..,7Vf-l,..)+ ^ v^a;/h 



(7iii) 



It can be noted that right-hand sides of the equations (7i) and (7ii) contain 
only boundary terms whereas the right-hand side of (7iii) contains many 
more than just boundary terms. 



3. Gauss's theorem and conservation laws in a discrete space 

Let a domain D of the four-dimensional discrete space be given by (see 
Fig. 2) 

D := {n e : iVf < n^" < N^, e {1, 2, 3, 4}} . (8) 
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n 



(Nl,N?) 



FIG. 2 A two-dimensional dicrete domain. 



The discrete boundary points of D are taken to be 

dD = di-D U di+D U d2-D U d2+D U da-D U da+D U a4-L' U 84+0, 
d^^D := {n e : = N^, <rf < N^, a ^ /x}, (9) 



We also denote the unit "normal" on the boundary dD by the follow- 
ing definition. 

{1 on da+D , 
-1 on . 



We assume that a tensor field j'^,"{n) is defined over D C N*^. (See equa- 
tion (34A) for the definition of a tensor field.) Now we are ready to state 
and prove formally the "discrete Gauss's theorem" [5] . 
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Theorem 3.1 (Discrete Gauss's): Let a tensor field j^"{n) be defined on 
DUdDcN"^. Then 

TVj 7V| iV| iV| (4) 

E E E E [^.j'/in)] =: E [A,j^-{n)] 

n^=Nlv?=Nln^=Nln-^=N^ DcN* 

(3) 

= E ^'■'^■■Wi^MW- 



Proof. Using the equation (7i) four times to the left-hand side of the 
equation (11) we obtain 



N'i JVl 

E E E 

n2=Ar2 n^=Nf n'^=Nf 



Ni Ni N* 



E 



+ .. + ..+ 



= E E E U'-(^2,-)-j?-(A^iS ■■)] + ■■ + ■■ + ■■ 

n2=Ar2 n3=Af3 n'^=Nf 
(3) (3) 

E E J.'-H 

di+D di-D 



(3) 

E j.^'W^iW 

di+Dudi-D 

(3) 



+ .. + ..+ 



□ 



We shall now make some comments on the preceding theorem. 

(i) An alternate theorem holds by replacing the left-hand side of the equa- 

(4) 

tion(ll)by E [A'^j^-in)]. 

DCN4 

(ii) Both forms of Gauss's theorem can be generalized to any finite-dimen- 
sional discrete space. 

(iii) In the finite difference representation of quantum mechanics [3], the 
four momentum operators are furnished by = — «A*. These are conse- 
quences of the relativistic representations. However, Gauss's theorem 3.1, 
which uses operators, is non-relativistic. (The relativistic Gauss's theo- 

(4) 

rem involving J2 [A*j'^- (n)] is not yet solved. See the comments at the end 
Den* 
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of this section.) The partial difference and difference-differential conservation 
equations (non-relativistic) are furnished by 

A^j^(n) = 0, (12A) 

A.j'rM+dtjt-M^O. (12B) 

Difference conservation equations lead to summation conservations. We shall 
presently state and prove a theorem about this topic. 

Theorem 3.2 (Conserved sums): Let a tensor field satisfy the 

partial difference conservation A^j^-^n) = in a domain D C given in 
the equation (8). Let furthermore the boundary conditions 

f::{n)vk{n)\Q^^Dvjda-D = (13) 

hold. Then the sum 

ATI Ar| Af| 

E E E \jt-{n\n\n\Nt)] 

n^=Nl n'^=Nl v?=Nf 
Nl Af| Af| 

- E E E \j'r{n\n\n\Nt)] (14) 

n^=Nl n2=Nl n^=Nl 
Nl Nl 7V| 

- E E E 3':-{n\n\n\n') 

n^=Nl 'n?=Nl n^=Nf 

for all satisfying Nf < < ■ 

Proof. Using Gauss's theorem 3.1 and the difference conservation equa- 
tion (12A), we conclude that 

(4) 

= E [A,j^:{n)] 

(3) (3) 

E £-(n)z/i(n) + .. + ..+ J^:\n)ui{n). 

di+Dudi-D dA+DVJdi-D 
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Assuming the boundary conditions (13), the above equation yields 

(3) 

di+DUdi-D 
(3) (3) 

di+D di-D 

Considering Gauss's theorem for a proper subset of D, the equation (14) 
follows. □ 



In the case of the difference conservation equation (12) being valid for 
the denumerably infinite domain N"^, we can derive conserved sums under 
suitable boundary conditions. Such boundary conditions (sufficient) are 

da-D := {n e : = 0, < n'' < oo, a 6} , 

da+D := {n e : = M, < n'' < oo, a 7^ 6} , (15A) 



lim 

M— >oo 



Ida-DUda+D 



0; 



da-D ■= {(n,t) G X i? : = 0, < < oo, a ^ 6} , 

da+D := {(n,t) e X : = M, < n'' < oo, a 7^ 6} , (15B) 



lim 

M— »oo 



j':in,t)iy,{n,t) 



\da-DUda+D 



0. 



Under boundary conditions (15A), the equation (14) yields the following 
totally conserved quantities (generalized charges!) 



oo oo oo 



^3:; = E E E 

ni=0 n2=0 n3=0 



(3) 

=: E J-.'-ln^^') 

nGN3 

(3) 

= Ei.'-(n'2). 

neN3 



(16A) 



In the case of the difference-differential conservation equations (12B) and 
the boundary conditions (15B), we can derive the totally conserved quantities 



(3) (3) 
nGN3 nGN3 



(16B) 



One may wonder why are we considering the non-relativistic Gauss's the- 
orem and the consequent conserved sums at all! In Section 5, we shall prove 
that the theorems of this section can he used tactfully to elicit relativistic 
conserved sums. 
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4. Relativistic covariance of pEirtial difference and difference- dif- 
ferential wave equations 



The relativistic invariance or covariance of our partial difference or difference- 
differential equations is quite delicate. The criterion used in this paper has 
been developed through many papers [6] published in the last three decades. 
Let us start with a very simple example of Poincare transformations and the 
invariance of the usual wave equation under that transformation. Consider 
the infinitesimal time-translation: 

a 4 ^4 4 V / 

There exists in the old frame a different event (x*) which has the same 
coordinates (x) in the new frame. The coordinates of (x*) from (17) are 
given by 

^ ^a^ ^#4 ^ ^4 _ ^4^ j#4 ^ ^4 _ ^^g) 

Consider the transformation rule for a scalar field (p{x) given by 
${x) = (j){x) , 

(f){x) = (l){x*) = (f){x^,x'^,x^,x'^ — e'^) . ^^^^ 

Let be a Taylor-expandible (or analytic) function. In that case we 
can express (19) by Lagrange's formula as 

0(a;) = [exp(— £^94)] (f){x) 

= [exp{—ie'^{—i)d4)] (j){x) 
= (P{x) - e*di<j>{x) + [{e^f] 

:- 



(20) 



^x^^ 

Moreover, let ^(x) satisfy the usual wave equation 

ri^'''d^dA{x) = Q. (21) 

In the new frame, 

r]'"'d^d^4){x) = r]^"'d^d^{[exp{-s^d4)](t){x)} 

= [exp(-£4a4)] [r''d^d,(P{x)] = . 

The above equation demonstrates in an unusual manner the relativistic in- 
variance of the wave equation under an infinitesimal time translation. We 
shall follow similar proofs in the sequel. 
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Let us now re-examine the very concept of relativistic invariance or co- 
variance. The relativistic covariancc docs not necessarily imply that the 
space and time coordinates must be treated on the same footing. Nor do the 
equations which treat space and time variables on the equal footing auto- 
matically imply the relativistic covariance. For example, let us consider the 
partial difference Klein-Gordon equation [7] in the lattice space-time as 

r/^'^A^A;0(n) - m^(t){n) = . (22) 

This equation treats discrete space and time variables on the same footing. 
However, the equation (22) is certainly not invariant under the continuous 
Poincare group TO {3, 1)! (Although there exists the lattice Lorentz group [8], 
a subgroup of the Lorentz group 0(3, 1), which leaves the lattice space-time 
invariant.) 

Consider another example, namely the Friedmann-Robertson- Walker cos- 
mological model of the universe. The corresponding metric and the "or- 
thonormal" tetrad are given by [9]: 



ds'^ = [R{t)]^[l + K5abx''x^]-^[5ijdx'dx^] - (dt) 
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ef.) = [R{t)]-^[l + KSabX-x']6^j^, (23) 

6(4) - ()(4) . 

The above metric and the corresponding tetrad do not treat the space and 
time variables on equal footing. (Although these are extracted as exact so- 
lutions of Einstein's general covariant equations.) However, if we consider a 
suitably parametrized motion curve given by a time-like geodesic, the appro- 
priate Lagrangian and the four-momentum are given by: 



L(..) = {m/2){[R{t)]-'[l + KSabX'^x']-%,x'x^]-ii^)} , 
PU) := ^^e'(.^^m[R{t)][l + K5abx''x']-'5jix\ 



(24) 
m > 0. 

Recalling that g^j^i,[x)x^x'' = —1 along a time-like geodesic, we obtain from 
(24) that 

r]^''^^''^P(,)Piu) = -m' . (25) 

Therefore, the special relativistic equation (25) holds among the tetrad com- 
ponents of the four-momentum locally. In fact, in every reasonable curved 
universe, with any admissible coordinate system, the special relativity holds 
locally whether or not the Poincare group IO{3, 1) is globally admitted as 
Killing motion. 
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In a flat space-time, the first quantization of the equation (25) leads to 
the relativistic operator equation: 



[r]'"'P^P, + m'l](l) = 0. (26) 

— * — * 

Here, and are vectors in the tensor-product [3] of the Hilbert spaces. 
The P^'s represent four self-adjoint, unbounded hnear operators and "1" 
stands for the identity operator. The equation (26) physically represents the 
quantum mechanics of a massive, spin-less, free particle. Let us explore the 
relativistic invariance of the operator equation (26). The finite and infinites- 
imal versions of a Poincare transformation in the classical level are given 
respectively by: 

(27) 

— _|_ £M _|_ £l^x'^ ^ 
Sfxu := Tl^aSl = -£un + 0(£2) . 

However, there exist in a deeper level, the quantum Poincare transfor- 
mations on quantum mechanical operators. In fact, there are two possible 
quantum Poincare transformations of operators which involve the same uni- 
tary mapping. The Heisenberg-type of Poincare transformations [10] under 
the infinitesimal version of (27) are furnished by: 

U{e) := exp {-ie'^P^ + {i / ^)e>''' {Q ^P, - Q,P^ + P,Q^ - P^Q,)} , (28iH) 

Q^.P.-PvQ^. = ^V^.^. (28iiH) 
P^ = U\e)P^U{e) = - £^P, + 0(£2) , (28iiiH) 

= U\e)Q^U{e) = + + e^g^ + Q{e^) , (28ivH) 
= 0. (28vH) 

Here, the dagger denotes the hermitian conjugation. 

In the Schroedinger-type of covariance [10], the abstract operators and 
the state vectors transform under the infinitesimal version of (27) as: 

P^. = P^.. Q^ = Q'', (28iS) 
$^U{e)$, (28iiS) 



5l$ := 0-0 = {-i[e^P^-{l/4)e>'- 

(Q^^P. - QuP^. + P.Q^. - PM] + O(£')}0 ■ 



(28inS) 
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Here, U (e) is the same as in (28iH). The vector 5l(P is called the Lie- variation 

of the vector (f). 

Note that in the two types of transformations, the expectation values of 
a polynomial operator F{P, Q) is given by: 

(|,|F(P,Q)||„)^, = {MF[U^PU,U^QU]\$a) 

= {l\F{P,Q)\l),. 

(Here, {(f)b\4'a) denotes the inner-product of the two vectors (pb and 0^.) 
Therefore, physical quantities transform exactly in the same manner under 
Heisenberg-type and Schroedinger-type of quantum covariance rules. We 
shall follow the Schroedinger-type of covariance in this paper. 

It is well known [11] that the operator P^Pfj,, which is one of the Casimir 
operators of the Poincarc group X0(3, 1), commute with all the generators 
P^ and {l/^){Q^Pu-QuP^i+PyQ^i-P,j,Qy) of the Poincare group. Therefore, 
we obtain from (28i-ivS), 

[P^P^ + I ] = [P'^P^ + I ] [0 + 5l$\ 

= [P-P, + mH][5L$] 

= -I \e^P^-{\l^)e^^{Q^.P.-Q.P^^P.Q^-PM\-x ^^^^ 
[P<^P<,+m2l]0+O(£2) 

= 0(£2) . 

Therefore, the operator equation (29) is an actual proof of the relativistic in- 
variance (up to the 2nd order terms) for the operator Klein- Gordon equation 
(26). 

We can generalize the Lie-variation (28ivH) for an arbitrary relativistic 

— * 

tensor (or spinor) operator - . The appropriate definition is given by: 

J- - = Sl$- := -i {e'^P^ - (l/4)£'^'^ ^^^^ 
[(Q^P. - Qt.Pu + PuQ^ - PM + 2S-^,,] }f + 0{e^) , 

where S'^^, = —S'^^ denotes the "spin operator". 

The vahies of the entries for S'^^, can be calculated exactly from the usual 
tensor and spinor transformation rules. For example, if we consider a vector 
field 0" with components (f)"', then 5"^^^ = rji^i^S'^ — rjn/sS". 

Now we shall state a necessary criterion for the relativistic invariance of 
an operator equation. 
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"Let an infinite- dimensional Hilhert space vector (f)" represent a particle 
with zero or non-zero spin. Let the corresponding Lie-variation ' he given 
by the equation (30). In the case of the operator equation for 0" being rela- 
tivistic invariant or covariant, the Lie-variation 6l4>" must satisfy the same 
operator equation as (f)" does {up to the 2nd order terms)." 

Let us apply the above criterion of covariance on different first-quantized 
equations arising out of different representations of quantum mechanics. 
Firstly, consider the usual Schrocdinger representation of quantum mechan- 
ics, namely, = —id^ and = x^. The operator equation (26) of the 
mass-shell constraint goes over into 

-U^{x) + rn^4){x) := -[r]'"'' d^d^4){x) - nv^^ix)] = . 

This is the usual Klein-Gordon equation. The Lie-variation under TO(3, 1) 
is given by 

5L(t>{x) = ^ [d'd^ - (1/4)5^'^ (x^a, - x,d^ + d^x^ - d^x,)] ct>{x) + 0(£2) . 
Therefore, 
— □(5i0(x) -\- m?5L(j){x) 

= -[e>'d^-{l/A)ei'-{x^d,-x,d^+d,x^-d^x;)\ [-□0(x)+m20(x)]+O(£2) 
= 0(6^) . 

Therefore, the relativistic invariance (up to the second order term) is assured. 

In the second example, let us consider the finite-difference representation 
[3] of the quantum mechanics by putting 

P^ = -iA*, (1/V2)(A^V^- Vr^A;, + 2V^I). 

Here, the index is not summed. The operator equation (26) yields 

77^'^A#A*0(n) - m20(n) = , 
n := {n^,n^, n^, n"^) G . 

This is the finite difference version of the Klein-Gordon equation. The cor- 
responding Lie- variation is given by: 

-(A^Vn^ - Vn^K + 2Vn^)A* + A*{A^^/n^^ - ^fnPA'^ + 2^) 
-A#(A^V^ - ^/nPK + 2V^)]0(n)} + 0(^2) . 

Here, the indices v are summed. It can be proved by direct computations 
that 

77'^'^A*A*[5^0(n)] - m''[5L(t>{n)] = 0(£^) . 
14 



Thus the partial difference equation (31A) is indeed invariant (up to the 
second order terms) under the continuous Poincare group! 

Now, as the third example, let us consider a mixed difference- differential 
representation of the quantum mechanics by choosing 



Pb = -iAf, = -id^ = idt , 



x'^ = t. 



Here, the index h is not summed. The operator equation (26) yields the 
difference-differential version of the Klein-Gordon equation: 



(5"^A# A*0(n, t) - {dtf(P{n, t) - mV(n, t) = , 
(n, i) e X K . 

Here 5"^ is the Kronecker delta. We have for the Lie-variation 

SL(j){n,t) = 

-{AbVnP - VnPA'f, + 2VnP)A* + Af (A"v^ - Vn^K + 2v^) 



(31B) 



(32B) 



-A*(A5^- v^A';, + 2v^)] (/)(n, t) 

+(1/ V2)£''^ [{AaV^ - V^A'Jdt - tAt] }0(n, t) . 

Here, the indices a, b are summed. We can prove by a long calculation that 

5«^A*A*[5^0(n, t)] - imS^Hn, t)] - m'[Sr.<f>{n, t)] = 0{e') . 

In other words, the difference-differential equation (3 IB) is indeed invariant 
(up to the second order terms) under the ten-parameter continuous group 
T0(3,l)! 

As the last example, consider the Schroedinger's difference-differential 
equation [3] for a free particle (with m > 0) as 

(2m)-M«^A#A#V'(n, t) + idti^{n, t) = . 

The Lie- variation ^^^/'(n, t) according to the equation (32B) does not satisfy 
the Schroedinger difference-differential equation up to the second order terms. 
The reason for this failure is the operator P4 = idt does not commute with 
three particular generators Q4P6 — QhP^ of the Poincare group. Therefore, 
wc obtain an alternate proof of the well-known fact that the Schroedinger 
equation is non-relativistic. Thus, our necessary criterion involving the Lie- 
variation can prove or else disprove the relativistic invariance or covariance 
(up to the second order terms) of a quantum mechanical system in any rep- 
resentation. 
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We need to define the exact transformation rules for a tensor or spinor 
field in partial difference and difference-differential representations. Recall 
from (27) that the finite Poincare transformation is given by 



(33) 



The exact transformation rules for the tensor and spinor fields under (33) 
are furnished by: 

$-{n) = (exp{-a''A* + (l/8\/2)(a;'^'^-a;'^'^) 

(A^Vn^ - v^A'^ + 2v^) A# 

-(A^V^- V^A'^ + 2Vri^)A* (34A) 
+A#(A^^A^ - Vni'A'^ + 

-A*(A,V^ - V^A', + 2^) + i2S-^^\ }) • 0-(n) , 



0--(n,i) = (exp{-a^A#-a49t-(l/8^/2)(u;«^-u;^") 



^AaVn" - v^A; + 2v^) Af 
-(AbVn^ - Vn^A; + 2Vn^)A* 

+A*(A„^A^ - aA^A;, + 2a/^) (34B) 
-A#(A,,V^ - V^A; + 2v^) + i2S'^(,.." 
-(i/2^) (a;«4 _ ^4a) [(A„Vri« - Vn"A'J9t 
-iA*+z5„4..]})-0-(n,t). 



Here, indices and a, h are summed and 5"^^ stands for the spin-operator 
for the field 0- (n) or 0- (n, i). 

An obvious question that arises is how the discrete variables n^.n"^, , 
transform from one inertial observer to another! We can recall that the inte- 
ger appears as the eigenvalue 2n^ + 1 of the operator (P'*)^ -|- {Q^Y- A 
particle in the corresponding eigenstate is located on a circle of radius 
-\/2r7/M-T in the p'^— x'^-th phase plane (see Fig. 1). However, for a relatively 
moving observer, according to the Schroedinger-type of covariance, the par- 

— * — * — * 

tide is in the state = U{e)ipni^. (See equations (28iS,iiS,iiiS).) But tp is 
not an eigenstate of the operator {P^'f + {Q^f = {P^f + {Qi'f. Therefore, 



16 



the particle appears in a fuzzy domain in the phase plane for the moving 

— * 

observer. But for the moving observer, there exists another state ip^^, :— 

7^ u(s)jfr,. = ^ suchjhat [(P'^y + = [(P'^y + (Q''y]^fn^ = 

{2n^ + l)i/j^f,. Therefore, t/j^^ is an eigenvector for the operator (P^y + {Q^y 
in the moving frame and the corresponding eigenvalue is 2n'^ + 1. Thus both 
observers have exactly similar discretized phase planes {see Fig. 1), although 
discrete circles in one frame do not transform into discrete circles in another 
frame! 

Now, we shall prove the relativistic invariance of the summation operation 

oo oo oo oo oo oo oo 

H H H and the sum-integral operation E E J2 I dt. Consider 

nl=0 n2=0 n^=0 n*=0 n^=0 n2=0 n^=0 M 

the state vector (f) representing a scalar particle in the abstract equations in 
(26) and (28i,ii,iiiS). Mathematically speaking, (p e Hi (S> H2 (S) ^ Hi, is 
the tensor product of Hilbert spaces [3,12]. Here, the Hilbert space Hn is 

— * 

acted upon by the operators Pn, etc. The square of the norm of is 
given by: 

11^11':= (010)- 
But from (28iH), (28iiS) we have 

\\$f^{$\U\e)U{e)$)^U\\\ 

where the dagger denotes the hermitian-conjugation. Thus, the norm 

is invariant under the unitary transformation induced by an infinitesimal 

Poincare transformation. 

Moreover, under a finite Poincare transformation which induces a unitary 

mapping, we can conclude that || || = \\(f)\\. 

In the finite difference representation of the scalar field [3] we have 



00 00 00 00 



|2._ V- V- i^(n\n\n\n^)\\ 



:= E E E E 

ni=0 n2=o n3=0 n^=0 



Moreover, in the difference-differential representation 

00 00 00 „ 
Ur--= E E E / \Hn\n'y,t)\'dt. 

ni=On2=On3=o jg 

00 00 00 00 

Therefore, by the invariance of ||0|| , both operations E E E E and 

ni=0 n2=0 n3=0 n*=0 



00 00 00 



J2 J2 J2 I dt must be relativistically invariant. 

ni=0 n2=0 n3=0 R 
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5. Variational formalism and conservation equations 



Let us consider a real- valued tensor field A^^^{n) in the difference representa- 
tion and A'^"[n,t) in the difference-differential representations. The action 
sum or the action sum-integral for such fields are defined respectively by: 



Nl ATI 



A{A--):^ E E E E 

n^=Nl v?=Nl ■n?=Nl n'^=Nf 

L (n; y" - ; y'f;-) 



(35A) 



A{A--):^ E E E 

ni=Afi n2=Nl n^=Nf ti 



=A-- (n,t), J/-- =A* , yX-=dtA-- 



(35B) 



The Euler-Lagrange equations under zero boundary variations for A''^^{n) 
or A'^-{n,t) are given by (see Appendix I): 



dL{.. 



0, 



(36A) 



dL{..) 



dL{.. 



dL{.. 



0. 



(36B) 



In the case of a complex- value tensor field (p^'^n) and 0'^- (n, t), the action 
sum and sum-integral are defined respectively by: 



Nl Nl Nl Nl 



A{<t>-,r-):= E 

n^=Nl n'^=Nf n''^=Nf n4=A 



w 



(37A) 



Nl 



Nl 



t2 



Air-,r'):= E E E / 

n^=Nl v?=Nl nS=Nf ti 



|p''-=,^.'^--(n,t),p''-=0'---(n,t) 



dt. 



(37B) 
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Here, the bar stands for the complex-conjugation. 
The corresponding Euler- Lagrange equations are: 



dL{. 



dp"-]. 



- A* 



dL{.. 



0, 



(38A) 



dL{. 



dp"-]. 

dL{..) 
dp"-].. 



A* 



dL{.. 



dL{..) 
dpl- 



0, 



-dt 



dL{..) 
dpi- 



0, 



(38A) 



(38B) 



dp^'-y 



dL{.. 



■b J I 



d, 



dL{.. 



{38B) 



In the derivation of above equations, the techniques of the complex- 
conjugate coordinates are used [13,14]. 

We shall now derive the partial difference and the difference-differential 
conservation equations for various fields (see Appendix II). 



a^t; + .. = o. 



(39Ai) 



^Vv" \{..,n"-l,..) 



. A*A'^- 



+|^-(A*^'^il(..,--M-W--)|.. 



A,T^ + dtT^ + .. = 0, 



At,n + dtT^ = 0, 



(39Aii) 

(39Bi) 
(39Bii) 



dL{..) 



AtA'^) 



d{L.. 



dVa. 



dy''- 6|(..,nf-l,..) 

I •(A*A-)l(..,„._i,..)-<5^L(..)|. 



(39Biii) 
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dL{.. 



(39Biv) 



TIM ■■=\h7 



dL{..) 



b\{..,nb-l,..) 



, (39Bv) 



(39Biv) 



4|- 



We can sum or sum-integrate the relativistic conservation equations over 
an appropriate domain of or x M. Using Gauss's theorem 3.1 and 
assuming boundary conditions similar to the equations (15A,B), we obtain 
the total conserved four-momentum: 



00(3) 00(3) 

-p, = E Ki^, n') + ..] = E 2) + ..] , 



(40A) 



n=0 



n=0 



00(3) 00(3) 



n=0 



n=0 



(40Bi) 



00(3) 00(3) 

H := -P4 = E [Tli^^ t) + -] = E [Tli^^ 0) + ..] , 



n=0 



n=0 



(40Bii) 



oo{3) 00 00 00 

E := E E E . 

n=0 ni=0n2=0n'^=0 

The total energy-momentum components are given by incomplete 
equations (40A), (40Bi,ii) for a general Lagrangian. However, in the fol- 
lowing paper, we shall derive exact equations for the Klein-Gordon, electro- 
magnetic, and Dirac fields. 

In case of a complex- valued field the difference and difference-differ- 
ential conservation for the charge- current vector is given by equations 
(A.11.9A,B) as: 

A^/H=0, (4lAi) 



j^{n) := ie 



2 

dL{..) 



dL{..) 



.^PM|(..,nM-l,..) 



+ (c.c.) , 



Ml- 



(41ii) 
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A,j\n,t)+dtf{n,t)^0, 



(41Bi) 



/(n,t) := ie 




dL{..) 

Qpa.. 



h\{..,n''-\,..) 



r-(n,t) 



+ (c.c.), 



(41Bii) 



/(n,t) := ie 



dL{..) 
dpl 



(n,t)-^-r-(n,t) 



(41Bm) 



Here, (c.c) stands for the complex-conjugation of the preceding terms and 
e = ^47r/137 is the charge parameter. 

Under appropriate boundary conditions (15A,B), we can derive the con- 
served total charge: 



. 00(3) 

n=0 



^P4|(n,l) ^/^4|(n,2) 



fn,ll 



• + (c.c.), (42A) 



Q 



le 



00(3) 



dL{.. 



r-(n,0) 



+ (c.c.) 



(42B) 



Note that equations (42A,B) are already exact. 

We shall now make some comments about the total conserved quantities. 

In the preceding section, we have deduced that the four-fold summation 
00 00 00 00 00 00 00 

S S S S s-nd the summation-integration S J dt are rela- 

ni=0 n2=0 n3=0 n^=0 ni=0 n2=0 n3=0 IR 

tivistically invariant operations. We have obtained the total four-momentum 
Pfj_ in the equations (40A), (40Bi,ii) by the four-fold summation and the 
summation-integration of the relativistic conservation equations (A.11.4A) 
and (A.II.4Bi,ii) respectively. Therefore, we claim that (the complete) P^^s 
are components of a relativisitic four-vector. Similarly, the total charge in 
equations (42A) or (42B) is relativistic invariant. 

The physical contents of the partial difference conservation (39Ai) and 
the difference-differential conservation (39Bi) are identical. However, the 
total four-momentum components in (40A) and (40Bi,ii) are quite dif- 
ferent inspite of the same notations! The reason for this distinction is that 
in equation (40A), the summation of the field is over a = const, "hy- 
perspace" . In covariant phase space, an = const, "hyper surf ace" implies 
that {tY -\- {pif — (2n^-|-l) = const. Therefore, in the t-p'^ phase plane, the 
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time coordinate t oscillates between the values — \/2n^ + 1 and -\/2n^ + 1. 
It is certainly not a t = const, "slice". However, in the case of equations 
(40Bi,ii), the usual t = const, "hypersurface" of the discrete phase space and 
continuous time is used. 



Appendix I: Euler-Lagrange equations 



We shall consider a two-dimensional lattice function / for the sake of sim- 
plicity. We firstly choose a two-dimensional discrete domain (see Fig. 2) 

D:={neN^: < Nf^K N^, 2 <N^-N^,Nf<n^<N^,Ae {1,2}}. (A.I.I) 

We shall follow the summation convention on capital indices which take val- 
ues from {1,2}. 

A real- valued lattice function / : C — > M is considered. The 



Lagrangian function L : D C N x 



X 



is assumed to be Taylor- 



expandible (or real analytic). (In all practical purposes, L is a polynomial 
function.) The action functional A is defined to be the double-sum 



A{f):= E E Hn'^y^yA)\y=fin),yA=A*f{n) , U {n\ U^) . (A.I.2) 



Now we define the variations of the lattice function / by 
Sf{n) := eh{n), 
S[A*f{n)] eA*h{n)^A*[Sf{n)]. 



(A.I.3) 



Here, e > is an arbitrary, small, positive number. Therefore, the variation 
of the action functional A is furnished by: 



^A{f) = E E {L{ri]y;yA)\y=f+sf,y^=A*f+s{A*f) 

n^Nl n'^=N'f 

-L{n;y;yA)\y=f^n),yA=A*fin)} 
-dL{..] 



= ^ E E 

'dL{ 



dy 



(A.I.4) 



h{n) + 



dy/ 



■A*h{n)] + 0{e^) 



\y=f{n),yA=A*f 
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The variational principle, which implies the stationary (or critical) values 
of the action functional, can be stated succinctly as 

The above equation yields from (A. 1. 4) that 



Afi iV| 

E E 



dL{..) 
dy 



h{n)+ 



J |(ni,n2 



dVA 



A*h{n) } = 0. (A.I.5) 



J |(ni,n2) 



(Here, we have simplified the notation by putting |(n^,n^) in place of \y 



Va 



A^^f {n^ , n'^) .) Using the definition (4iii), opening up the 



double-sum in (A.I.5), and rearranging terms, we obtain (after a very long 
calculation) 



E E 



1 


ldL{..)] 


-A* 

^A 


raL(..)i 


|(ni,n2)} 




dy 


|(ni,n2) 







h{n\n') 



(A.I.6) 



+ (Boundary terms) = . 



Here, the boundary terms are given by: 



(Boundary terms) := 

[n! ""^^^ dLj..) 
V 2 dy2 \{Nl+k,Nl) 



h{Nl + k, - 1) 



+ 1 



2 k=0 \{N^+l,N^+k) 

NjTl dLj..) 

2 k=0 \{Nl+k,Nl) 



h{N^ + l,N^ + k) 
h{Nl + k, Nl + 1) 




Nl-N^ 

E 



dL{..) 



+ E 



2 k=o miNi+k) 

dL{..) 



dy \(Ni+j,Ni) 



h{Nl-l,N^ + k) 
dL{..)\ 



A? 



' N? + 1 dL{..) 

2 dy2 \^Nl+j,N^+i) ^ 
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V J\.. 
h{Nl + j,N^) 



+ E 



dL{..) 

9y \{NiNi+j) 



|(JVi,iV2+j)J 



^^^-^^-^ I dLi..) 



+ E 

i=i 



2 9yi |(^i_i^^2^^-) 



\{Ni+j,Ni) 



A* 



dy 



1 |. 



\{Nl+j,Nl) 



+ 



'iV| dL{..) 



2 ^?/2 |(iVi+j,iV|) 

+ E 



■h{Nl+j,N^) 



dL{..) 



iNl + 1 dL{.. 



+ 



j=i I 9y \{Nl,N^+j) V 2 dyi I (^1+1^^2 



(9L{.. 



dy2 



I- 



|(iVi,iV2+i)J 

1 



(^y i(ivi,iv2) \/2 

aL(..) 



\/nI + 1 



+ 



^^2 |(JVi,iV2+l) 

Kivi.iv?) \/2 

^^(■■) 
9y2 |(Ari,Ar2+i) 



+ 



dL{..) 



1 

+ ^ 



9y \{N^,Ni) a/2 
9L(..) 



+ 



dy2 |(ivi,Ar|-i)_ 

(^y i(ivi,jv|) \/2 

9L(..) 



\{Nl,N^-l) 



dyi \{Ni^i,Ni) 
KNlNl) 

T dL{..) 

9yi |(JVi-i,Ar2) 

■h{N^,Ni) 
- dL{..) 



dyi |(ivi-i,jv|) 
h{N^,N^) 



\/nI + 1 



dL{..) 



9yi \{N^+i,Ni) 

h{NlN^). 
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(A.L7) 



The imposed plus natural boundary conditions [13] imply that the bound- 
ary terms in (A. 1. 7) add up to zero. There exist many possible boundary 
conditions. We shall adopt the simplest of all, namely 

(A.I.8) 

5/ (n^n^) I Boundary = 0. 

Note that the discrete boundary points in equation (A. 1. 7) are more in 
number than those in equation (9). We have to prove now the analogue of 
the Dubois- Reymond lemma [14]. 

Lemma: Let the double sum 

Y: E 9{n\n')h{n\n')^0 

for a function g and an arbitrary function h over the domain D C defined 
in equation (A.I.I). Then gin^^n^) =0 in D. 

Proof. Choose h{n\n^) := ^^i'^^' for some {m\m'^) e D. (The 5^ 
denotes the Kronecker delta.) Then 

E E 9{n\n')h{n'y) 

Afi-l iV|-l 

= E E g{n\n')S!::s:g ^g{m\m')^0. 

ni=Afi+l n'^=Nf+l 

Since the choice of (m^,m^) & D is arbitrary, it follows that g{n^,n'^) = 
for all (n\n2) e □ 

At this stage we have furnished essentially the proof of the following 
generalization of the Euler-Lagrange theorem. 

Theorem (A.I.I): Let a function / : D C — > M, where D is defined in 
(A.I.I). Let an action functional A{f) be defined by (A. 1. 2). The stationary 
values of A, under the boundary variation 5f{n)\__ — 0, are given by the 
solutions of the partial difference equation: 



\y=f{n),yA=A*f{n) 



A* 



dL{. 



= . (A.I.9) 



|j/=/(n),J/A=A*/(n) 



(Here, the index A is summed over {1,2}.) 
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Appendix II: Peirtial difference and difference-differential conser- 
vation equations 



Consider the discrete domain D gN"^ given by equation (A.I.I). Let y4"' (n) 
be a real- valued (r-l-s)-tli order tensor field over D C N^. (See equations 
(34A,B).) The Lagrangian L(y"";y"") is a real-valued, analytic, and rela- 
tivistic invariant function over a domain of the continuum R'^''^" x R'^'^^'^'' . 
The action functional A{A°'^^) is given by the four-fold sum (compare equa- 
tion (A.I.2)) 



^1 



7V| 



7V| 



E E E E 

n^=Nl n^=Nl n^=N'f n*=N^ 



(A.II.l) 



In the usual case, L is a polynomial function of the 4^^''+*)+^ real variables. 
The Taylor expansion of L is given by 



^dL{.. 



L{y--,y^-) + 



Qya.. 



1 

+2 



Qya..QyP. 



Ay- + f^A,- 



(Ay--)(A/-) + 



+ 



dy^-dy^- 



(A<-)(Ay^-) 



dy^'-dyS 
+ ... . 



(Ay"-)(Ay^- 



(A.II.2) 



Here we have followed the summation convention on every repeated Greek 
index. Moreover, for a second- degree polynomial function L, the additional 
terms denoted by . . . exactly vanish. 

Now, we investigate the partial difference operations on L{..). Using equa- 
tions (4iii), (A.II.2), (5iv), (6i,ii,iii,iv) and not summing the index /x, we 
derive that 



\/2 A* [L(y"- ; <-)|.ya..=^a..(„)_j^a..=A*A° 

= {^/nF+l) ■ L[A''■■{..,n^' + l, ..)] 2-^l'\y/^iF+l A''-{..,n''+l, ..,n^ +1, ..) 

-^/nyA'^■■{.., nf" + 1, .., - 1, ..)] } 
-(x/n^) • L{/l"-(..,n'^ - 1, ..); 2-^l'^[^Jn'' + lA''■■{..,n^' - 1, ..,n^ + 1, ..) 

-v^A"- (.., ni" - 1, .., - 1, ..)] } 
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A^-{n) - A;^A"- ; A*A"-(n) - A'^A#A° 



+ 



"aL(.. 



I- 



1 a2L(.o_ 
■^1 



^ • (A,A*A-) (A,A*A/^-) 
2 dy^ -dy^ 



+ 



-^/^{L{..)+ 



A/ 4a.. I ^ ^^^(-O A' /1"-VA' _ 



+ 



-A' A*^"- + --^-^^ -(A' A*^"-)(A' A*^'^") - .. 1 



= ^[L(..)]A*(1) + 



dL{.. 



I- 



1 a^Lf 
+ - ^ 



(A^A*^"-) + (A;A*>1"--) 
v/;^^(A^^--)(A^A^") 



2 dy'^-dyf^- , 

-v^(a;,^--)(a;^/5-) 



^1 



2 (ay--) (ay^-) 
-v^(a;^a*a--)(a;a*a'^ 



V^^^(A^A*A"-)(A^A*^'5-) 
+ .. 



dL{..) 



JL(..)] A*(l) + . [A*A-- - ^-■(n)A^(l)] 

• [A*AJA-- - (A;rA"-)A*(l)] 

-yl'^-(n)A#A"- + ^"■■(n)>l^-(n)A#(l)] 



+ - 



1 a2L(..) 



2 dy^ -dyt 



■ [A*( A;r^-- . A? A^-) - {I^tA--){/^t^tA^-) 



ii^tA^-) (A#A*^--) + (A#A--) (A#A^-)A*(1)] + .. . 



(A.II.3) 
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Dividing this equation by -\/2, bringing the left-hand side term to the right- 
by the Euler-Lagrange equa- 



dL(..) 



hand side, and equating ^ 
tion (36A), we finally obtain: 



A* 

V 



dL(..) 



,A#^"- + ^^i:l A*A*^"--A*fL(..)i 



-A°-(n)A#A/^-- - (A*>1°-) {AP-{n))] 
1 ffT( ^ 



+[A*(1 



9L(..) 



Qya.. 



I- 



+ 



1 a^Lf..! 



2 dy^ -dyty 



:a*a"-) (a*a^- 



+ 



0. 



(A.II.4A) 



Here, the indices a, /3, i/, o" are summed. We claim that the above equations 
involving the relativistic operator A* = iP^ constitute relativistic conserva- 
tion equations in the finite difference form. The corresponding relativisitic 
difference-differential conservation equations are 



+dt 



' dL{..] 



|..J 



Qya.. 
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4^^r#t •[Af(ArA--.A*A^-) 

-(A* A--) (A*Af A^^-) - (A*A'^-) (A*A*A--)] 



-(A* A--) (A*5,A«-) - (a.A'^-) (A*A*A--)] 



+ 



dyf |.. " ^yf- |.. 



2 dy^-dyS I 
dy^-dy^-^ 



Ai^-^ ■ (dtA^-) (dtA^-)] + . . . ^ , (A.II.4Bi) 



2 dyf-dy. 



4 |. 



. \dyt 



+ dt\ 



'dL{. 



I- 



I- 



6 I- 



atA*A"- 



+W| •(^*)'^""-W-)i..]} + ---o. 



(A.II.4Bii) 



The above equations containing operators A* = iPj and dt = iP^ on the 
same footing are relativistic. It is extremely hard to put (A.II.4A) and 
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(A.II.4Bi,ii) into the relativistic conservation equations: 



A*T'' = . 



A*Tl + dtT^ = , A*Tl + dtT^ = . 

However, using the equation (5v), the relativistic conservation equation 
(A.II.4A) can be cast into the form: 



U '^Vv |(..,n--l,..) 



+ ^ - • (A^^ )\{..,n--l,..) 



Qya. 



1 + ..=: A4r»] + .. = 0. 



(A.II.5A) 



Here, neither T^{n) are relativistic tensor components, nor A^^ is a rela- 
tivistic difference-operator. However, the combination A,^T^ -|- .. are 
components of a relativistic covariant vector! 

Similarly, from the relativistic difference-differential conservation equa- 
tions (A.II.4i,ii) we derive 



dL{..) 



\{..,n>>-l,..) 



(l/x/2)Ab|V;^ 

+ ^-(Af^"-)|(..,n.-l,..)-W-) 

• AfA-j + .. =: A,[T„^] + d^[Tt] + .. = , 



(A.H.SBi) 



(l/v^)A 



dL{.. 



^?/r|(..,n^-i,., 

) 



(A.H.SBii) 



i..U: A6r4^ + atr4^ = o. 
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We can generalize conservation equations (A.II.5A) and (A.II.5Bi,ii) for 
a complex- valued tensor or spinor field and 0"(n, to the following 
equations: 



(1/V2)A,< 
dL{..) 



dL{..] 



\{-,n--l,..) 



• (A*0-)|(..,„._i,..) + (c.c.) - 5;[L(..)l 

..=: a.t; + .. = o, 



-/fA;[L(..)],. 



(A.II.6A) 



:i/x/2)A 



bSn 



dL{..) 



•(A„V)|(..,n^-l,..) + (C-C-)-<^^[i^(-)]| 
+9* 



% l. 



dpi \.. ' ^1 
= A67;^ + atr„4 + .. = o, 



+ .. 



(A.IL6Bi) 



:i/V2)A,|Vn'' 
, dL{..) 



dL{..) 



dp. 



%|(..,n''-l,..) 
• dt(t>\{..,n^-l,..) + (c.c.) 



(A.II.GBii) 



=: A,T4^ + airf = 0. 



Here, (c.c.) indicates the complex-conjugation of the preceding terms. 

Now, we shall investigate the gauge invariance of the Lagrangian for a 
complex- valued tensor or spinor field 0"' (n) and 0"' (n, t) and the cor- 
responding difference and difference-differential conservation equations. A 
global, infinitesimal gauge transformation is characterized by: 
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<p"-{n) = [exp(i£)]0"-(n) = 0"-(n) + (i£)0"-(n) + 0(£^), (A.II.7A) 
0"-(n,t) = [exp(i£)](/."-(n,t) = (/)"-(n,t) + (i£)0"-(n, t) + 0(e2). (A.II.7B) 



The invariance of the Lagrangian function under (A.II.7A) imphes that 



^(^■■'^" ;^M'^M)|p..=J..(n),p■■=J•■(n),p^=A*?■,p^=A*?■■ 

^ ( P' ■ , P' ■ ; , P;^ ) I p. . =0. . (n) , 75- ■ =0 ■ (n) , = A* ^- • , = A* 0- ■ 



(A.II.8A) 



its) 



dLL.) dU..) dU..) dLi..) 



dp- 



dp- 



Dividing by £ > 0, taking the hmit £ ^ 0+, and equating 



A* 



dH..) 



we obtain that 



A* 



' dL{..] 
. 5P^- 



|..J 



+ (c.c.) . (A.II.9A) 



Since iA* = Pp represent the relativistic four-momentum operators, the 
equation (A.II.9A) incorporates the relativistic partial difference equation 
for the charge- current conservation. The difference-differential version of the 
relativistic charge- current conservation is furnished by: 



Af 



' dL{..] 
. dpt 



r-(n,i) + 



dpt- 



r-{n,t) 



I- 



• A*r-{n,t) + 1^ • d,r-{n,t)] + (c.c.) = . 



(A.II.9B) 



dp^ 



b \. 



dpf 
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FIG. 1 One discrete phase plane. 



(N^Ni) 



(Nl,N?) 



FIG. 2 A two-dimensional dicrete domain. 



